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Abstract. The coupled numerical weather model WRF-SPA 
(Weather Research and Forecasting model and Soil-Plant- 
Atmosphere model) has been used to investigate a 3 yr time 
series of observed atmospheric CO 2 concentrations from a 
tall tower in Scotland, UK. Ecosystem-specific tracers of 
net CO 2 uptake and net CO 2 release were used to investi- 
gate the contributions to the tower signal of key land cov- 
ers within its footprint, and how contributions varied at sea- 
sonal and interannual timescales. In addition, WRF-SPA sim- 
ulated atmospheric CO 2 concentrations were compared with 
two coarse global inversion models, CarbonTrackerEurope 
and the National Oceanic and Atmospheric Administration’s 
CarbonTracker (CTE-CT). WRF-SPA realistically modelled 
both seasonal (except post harvest) and daily cycles seen 
in observed atmospheric CO 2 at the tall tower (R 2 — 0.67, 
rmse = 3.5ppm, bias = 0.58 ppm). Atmospheric CO 2 con- 
centrations from the tall tower were well simulated by CTE- 
CT, but the inverse model showed a poorer representation 
of diurnal variation and simulated a larger bias from ob- 
servations (up to 1.9 ppm) at seasonal timescales, compared 
to the forward modelling of WRF-SPA. However, we have 
highlighted a consistent post-harvest increase in the sea- 
sonal bias between WRF-SPA and observations. Ecosystem- 
specific tracers of CO 2 exchange indicate that the increased 
bias is potentially due to the representation of agricultural 
processes within SPA and/or biases in land cover maps. The 
ecosystem-specific tracers also indicate that the majority of 
seasonal variation in CO 2 uptake for Scotland’s dominant 
ecosystems (forests, cropland and managed grassland) is de- 
tectable in observations within the footprint of the tall tower; 
however, the amount of variation explained varies between 


years. The between years variation in detectability of Scot- 
land’s ecosystems is potentially due to seasonal and inter- 
annual variation in the simulated prevailing wind direction. 
This result highlights the importance of accurately represent- 
ing atmospheric transport used within atmospheric inversion 
models used to estimate terrestrial source/sink distribution 
and magnitude. 


1 Introduction 

The global climate is changing and these changes are driven 
by human activities, in particular by anthropogenic emis- 
sions of CO 2 (IPCC, 2007). The terrestrial biosphere cur- 
rently absorbs a significant fraction of anthropogenic emis- 
sions of CO 2 (Canadell et al., 2007). However, terrestrial 
ecosystems are highly complex and dynamic, creating a net 
land-atmosphere surface exchange that can be either a source 
or sink of CO 2 . Furthermore, the magnitude of sources and 
sinks vary both spatially and temporally, resulting in signif- 
icant seasonal, interannual and spatial variation. The mag- 
nitude of net ecosystem exchange of CO 2 (NEE) is signifi- 
cantly impacted by changes in weather, climate and human 
management, adding further complexity to ecosystem pro- 
cesses (IPCC, 2007). Higher spatial and temporal resolution 
observations over multi-annual periods are required to detect 
fine scale ecosystem heterogeneity and ecosystem response 
to drivers. A critical objective is an improved understanding 
of the information contained within observations of atmo- 
spheric greenhouse gas concentrations. 


Published by Copernicus Publications on behalf of the European Geosciences Union. 


736 


T. L. Smallman et al.: Investigating atmospheric CO 2 observations using WRF-SPA 


Both forward running and atmospheric inversion mod- 
els have been used in conjunction with observations of at- 
mospheric CO 2 concentrations to investigate regional scale 
exchange of CO 2 . Forward running models, such as WRF 
(Weather Research and Forecasting model) or RAMS (Re- 
gional Atmospheric Modeling System), have been used to 
investigate a wide range of topics, including the impact of 
surface processes on atmospheric transport (e.g. Steeneveld 
et al., 2011) and how ecosystems contribute to observations 
made at the regional scale (e.g. Tolk et al., 2009). Inverse at- 
mospheric models, which infer surface fluxes from measure- 
ments of atmospheric CO 2 concentrations (e.g. made at a tall 
tower), have been used to constrain the terrestrial carbon bal- 
ance at global, continental and regional scales (Gurney et al., 
2002; Peters et al., 2010; Lauvaux et al., 2012). Inverse mod- 
els are able to detect large-scale, large magnitude interannual 
variations in CO 2 exchange. For example, the Europe-wide 
heat wave in 2003 has been linked to a large-scale reduc- 
tion in carbon sequestration across Europe (Ciais et al., 2005; 
Peters et al., 2010). However, there remains uncertainty over 
the ability of regional scale inversions to successfully quan- 
tify small magnitude changes in surface fluxes. Peters et al. 
(2010) suggested two contrasting hypotheses to explain in- 
creases in regional estimates of sequestration in the years 
following the European 2003 heat wave. First, increased se- 
questration may have been due to interannual variation in 
plant phenology. Second, changes to atmospheric transport 
due to interannual variation in weather (e.g. due to turbulent 
exchange) may have significantly altered the footprint of tall 
tower observations. 

Observations of atmospheric CO 2 concentrations made at 
tall towers contain seasonal and interannual phenological in- 
formation about ecosystems near the tower (Miles et al., 
2012). However, the correlation between surface NEE and 
observed profiles of atmospheric CO 2 concentration declines 
with increasing distance from the observation tower (Gerbig 
et al., 2009; Miles et al., 2012). A reduction in correlation 
between tall tower observations and ecosystem activity is 
consistent with signal dilution due to atmospheric transport 
(Gerbig et al., 2009; Miles et al., 2012). The dominant in- 
fluence on observations of atmospheric CO 2 concentrations 
are from the near field, although the total footprint can cover 
a large area. For example, an area of ~500kmx 700 km 
has been simulated to contribute up to ~50% of the ob- 
served signal at the Cabauw tall tower in the Netherlands, 
while the land surface at the edge of this area contributes 
~ 10 time less than the land surface directly beneath the 
Cabauw tower (Vermeulen et al., 2011). Similarly, Gerbig 
et al. (2009) investigated the Harvard Forest tower, USA, 
estimating that the land surface fluxes within 20-60 km 
of the tower contributed a similar amount to observations 
as all other areas within their simulated domain combined 
(5000 km x 5000km). Furthermore, atmospheric inversions 
are unable to attribute variations and/or anomalies in atmo- 
spheric CO 2 concentrations to a specific ecosystem process 



Fig. 1 . Land classification map used covering the spatial extent 
of the model domain. The left panel is the parent domain at 
18x18 km, right panel is nested domain at 6 x 6 km. The star indi- 
cates the location of tall tower Angus. The map used in WRF-SPA is 
a modified MODIS land cover map provided with the WRF model. 
The fractions of each land cover within the nested domain are crop 
36 %, evergreen forest 1 %, mixed forest 42 %, grassland 2 %, man- 
aged grassland 13 %, upland 3 %, urban 3 %. 


(e.g. respiration or photosynthesis). A forward model has the 
advantage that processes and source area contributing to the 
atmospheric CO 2 signal can be directly investigated. 

We used a mesoscale model, WRF-SPA (Weather Re- 
search and Forecasting model coupled with Soil-Plant- 
Atmosphere model), operating in forward mode to simulate 
a 3 yr period between 2006 and 2008 over northern Britain 
(Fig. 1), which includes Scotland, the geographical focus of 
this paper. Simulated atmospheric CO 2 concentrations from 
WRF-SPA and a global atmospheric inversion model were 
compared to observations made at tall tower Angus (TTA) 
on the east coast of Scotland. TTA is currently Scotland’s 
only tall tower equipped for measurement of atmospheric 
CO 2 concentrations. WRF-SPA provides a means to upscale 
land surface exchanges, such as photosynthesis and respira- 
tion, using atmospheric CO 2 tracers, to observations made 
mostly within the planetary boundary layer (PBL) of region- 
ally integrated CO 2 concentrations, i.e. those made at TTA. 

The overall aim of this paper is to use ecosystem-specific 
CO 2 tracers of net uptake and net release of CO 2 to improve 
our understanding of how different ecosystems contribute to 
observations of atmospheric CO 2 concentrations and relate 
these contributions to surface processes (Tolk et al., 2009). 
A unique aspect of this study is the 3 yr observation data set 
used to consider the detection of both seasonal and interan- 
nual variation. 

Here we address the following questions: 

i. Does WRF-SPA more accurately simulated observed 
atmospheric CO 2 concentrations compared to a coarse 
resolution global atmospheric inversion model? 


Biogeosciences, 11, 735-747, 2014 


www.biogeosciences.net/ll/735/2014/ 



T. L. Smallman et al.: Investigating atmospheric CO 2 observations using WRF-SPA 


737 


ii. Can ecosystem-specific CO 2 tracers be used to inform 
on which ecosystem processes and land covers are re- 
sponsible for observed variations in atmospheric CO 2 
concentrations? 

iii. Can observations made at TTA detect variation in 
ecosystem carbon uptake, for ecosystems within 
the footprint of TTA, at seasonal and interannual 
timescales? 

2 Model description: WRF-SPA 

WRF-SPA (Smallman et al., 2013) is a coupling between the 
high resolution non-hydrostatic mesoscale model Weather 
Research and Forecasting (WRF) and the mechanistic land 
surface model (LSM) Soil-Plant-Atmosphere (SPA). SPA is 
fully integrated into the WRF model framework, where WRF 
simulates meteorological fields and atmospheric transport of 
the CO 2 fields. SPA in return provides WRF with surface 
temperature, roughness length, albedo, and exchanges of en- 
ergy, heat, water and CO 2 . A brief description of the WRF 
and SPA models are given below. 

2.1 WRF 

The Weather Research & Forecasting model (WRFv3.2) 
(http://www.mmm.ucar.edu/wrf/users/, accessed 19 Octo- 
ber 2009, 15:00 UTC) is a state-of-the-art non-hydrostatic 
mesoscale meteorological community model (Skamarock 
et al., 2008). WRF provides a highly adaptable model frame- 
work into which SPA has been integrated. In addition to sim- 
ulating meteorology, WRF is responsible for simulating at- 
mospheric transport of CO 2 tracers released by SPA and orig- 
inating from the anthropogenic and oceanic flux maps as well 
as lateral boundary and initial conditions. WRFv3.2 includes 
a number of land surface maps that include vegetation type 
and soil classification used by SPA and also orography, which 
impacts simulation of air flow within the model (Mesoscale 
and Microscale Meteorology Division, 201 1). 

2.2 SPA 

The SPA model is a high vertical resolution mechanistic ter- 
restrial ecosystem model (up to 10 canopy layers and 20 soil 
layers). SPA provides WRF with surface fluxes of heat, wa- 
ter and CO 2 exchange in response to meteorological drivers 
through a close coupling of its hydrological and carbon cy- 
cles, based on eco-physiological principles (Williams et al., 
1996). Detailed descriptions of the major SPA developments 
can be found in Williams et al. (1996, 1998, 2001, 2005), Sus 
et al. (2010) and Smallman et al. (2013). 

WRF provides SPA with meteorological drivers, includ- 
ing air temperature, precipitation, vapour pressure deficit 
(VPD), wind speed, friction velocity, atmospheric CO 2 mix- 
ing ratios, air pressure, and short- and long-wave incoming 


radiation. SPA currently has parameters for 8 vegetation 
types (evergreen forest, deciduous forest, mixed forest, 
arable cropland, managed grassland, grassland, upland and 
urban) suitable for UK application and 13 soil types im- 
pacting soil hydrology. Vegetation cover is specified by the 
MODIS land cover map provided with WRFv3.2, while soil 
classifications are from the default WRF soil cover maps 
(Mesoscale and Microscale Meteorology Division, 2011). 

The Farquhar model of photosynthesis (Farquhar and von 
Caemmerer, 1982), the Penman-Monteith model of leaf tran- 
spiration (Jones, 1992) and the leaf energy balance are cou- 
pled via a mechanistic model of stomatal conductance. Stom- 
atal conductance is modelled by linking atmospheric demand 
for water and available water supply from the soil through 
plant hydraulics (Williams et al., 1996, 2001). SPA max- 
imises carbon assimilation per unit nitrogen within a min- 
imum leaf water potential constraint to prevent cavitation 
(Williams et al., 1996). SPA uses detailed multi-layer param- 
eterisations of canopy processes, including radiative transfer 
(Williams et al., 1998), above and within canopy momentum 
decay and leaf level boundary layer conductance (Smallman 
et al., 2013). 

Plant phenology is described by a box carbon model to 
simulate the main ecosystem carbon (C) pools (Williams 
et al., 2005; Sus et al., 2010). C pools are foliage, structural 
wood carbon, fine roots, labile, soil organic matter (SOM) 
and surface litter. Crops have two additional C pools: storage 
organ C (i.e. harvestable C) and dead foliar C (still stand- 
ing). The C pools within WRF-SPA are “spun-up”, as de- 
scribed in Smallman et al. (2013), using meteorology, which 
is broadly representative of the median meteorological con- 
ditions in Scotland. The carbon model provides a direct cou- 
pling between the plant carbon cycle and plant phenology, 
specifically foliar and fine root C. Foliar C determines the 
leaf area index (LAI) while fine root C impacts water uptake 
potential. 

2.3 Initial and lateral boundary conditions for CO 2 

Initial conditions (IC) and lateral boundary conditions (LBC) 
for atmospheric CO 2 (except 2008) are from CarbonTrack- 
erEurope (CTE, Peters et al., 2010) providing 1° x 1° res- 
olution fields at 3-hourly update intervals. Optimised CO 2 
fields were not available for 2008 from CTE; instead, CO 2 
fields for 2008 are from CarbonTracker (Peters et al., 2007). 
These are also available at 3-hourly update intervals, but at a 
coarse 3° x 2° resolution. IC and LBC for atmospheric CO 2 
were linearly interpolated to the WRF-SPA domain. 

Global flux maps of anthropogenic CO 2 emissions and 
ocean absorption were used to provide non-biospheric sur- 
face CO 2 exchange. The global flux maps were also from 
CTE at 1° x 1° resolution with 3-hourly update intervals. 
Fluxes were interpolated using 4 point weighted mean based 
on latitude and longitude coordinates. Biospheric fluxes of 
CO 2 are simulated by SPA. All surface CO 2 fluxes were 
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calculated as rates which were added to the lowest model 
atmospheric layer in each time step. 

2.4 Atmospheric CO 2 tracers 

WRF-SPA has been modified with the addition of several at- 
mospheric CO2 tracer pools (Table 1). CO2 tracer transport is 
simulated within the model domain concurrently with meteo- 
rological variables; feedback on atmospheric radiative trans- 
fer due to variable CO2 is neglected as its impact is expected 
to be small. 

We compared simulated atmospheric CO2 with an active 
biosphere versus simulated without a biosphere (“forcings 
only”). “Forcings only” CO2 tracer contains IC, LBC, an- 
thropogenic emissions and ocean sequestration of CO2 (i.e. 
CO2 exchange not calculated by WRF-SPA) but no exchange 
with the SPA simulated biosphere. Comparison between the 
total atmospheric CO2 concentration and “forcings only” 
CO2 concentration allows for isolation of the impact due 
to inclusion of the simulated biosphere. Furthermore, CTE 
(2006-2007) and CT (2008) atmospheric CO2 concentra- 
tions are compared to TTA observations to assess how well 
CTE-CT simulates atmospheric CO2 at a tower not included 
in the inversion model. 

2.4.1 Ecosystem-specific tracers 

The ecosystem-specific tracers of net uptake and net release 
of CO2 are used to investigate the information content on 
these processes contained within the total atmospheric CO2 
concentrations simulated at TTA. We investigate how rep- 
resentative observations at TTA are of the underlying sur- 
face fluxes of CO2. Note that LBCs for the outer domain 
have been set with zero inflow and zero-gradient outflow for 
ecosystem-specific net uptake and net release CO2 tracers. 
Zero gradient inflow/outflow allow tracers to easily leave the 
domain and prevent artificial influx to the CO2 tracer fields. 

The land surface can be either a net source or net sink 
of atmospheric CO2, varying both spatially and temporally. 
Whether the land surface is a net sink or source of CO2 is de- 
termined by the net result of photosynthetic and respiratory 
processes. Atmospheric CO2 concentrations represent a spa- 
tial and temporal integration of the net flux of CO2 between 
the land surface and the lower atmosphere. Therefore, when 
the simulated surface flux of CO2 represents a net removal 
of CO2 from the atmosphere, an ecosystem-specific “net up- 
take CO2 tracer” is released into the simulated atmosphere at 
the same rate as the “surface net CO2 uptake flux” (i.e. rate 
of NEE). Correspondingly, when the surface net CO2 uptake 
flux represents a net addition of CO2 to the atmosphere, an 
ecosystem-specific “net release CO2 tracer” is released. 

The net uptake CO2 tracers are considered to be non- 
interacting/non-interactive, while the net release CO2 tracers 
are interacting/interactive. The net uptake CO2 tracers are 
non-interacting as they represent a removal of atmospheric 


CO2 and as such cannot interact with the land surface. After 
their emission from the surface, the net uptake CO2 tracers 
are transported through the model atmosphere. Conversely, 
net release CO2 tracers represent an addition of a physical 
mass of CO2 to the atmosphere via respiration, which can 
therefore be subsequently removed from the atmosphere by 
photosynthesis after its initial release. Allowing the removal 
of a net release CO2 tracer prevents a respiratory signal from 
being simulated at TTA, which in reality does not reach TTA 
due to being consumed en route in a physically consistent 
manor. 

Net release CO2 tracers are removed from the atmosphere 
if they are present in the lowest model atmospheric level and 
the land surface below represents a net removal of CO2. If 
there are multiple ecosystem-specific net release CO2 tracers 
present in the same model atmosphere grid box, then removal 
is determined by the relative fraction of each ecosystem- 
specific tracer. For example, to determine the removal of crop 
net release CO2 tracer. 


Cr ° P t CropT t fforest+ t %ass+ t Mother ’ 

where y f t cmv rm is the fraction of surface CO2 flux (i.e. 
NEE) to be applied to the crop net release CO2 tracer. f r crop 
is the crop net release (arrow indicating direction and t indi- 
cating it as a tracer) CO2 tracer concentration, similarly for 
forest, managed grassland and “other” land cover types. 

2.4.2 Investigating representativeness and detection of 
seasonal and interannual variation using CO 2 
tracers 

Ecosystem-specific CO2 tracers are used to infer representa- 
tiveness of the simulated atmospheric CO2 concentrations at 
TTA of surface CO2 flux and to investigate how much sea- 
sonal and interannual information is contained within these 
atmospheric CO2 concentrations. We assumed that the sim- 
ulated atmospheric CO2 tracers are driven by the simulated 
surface CO2 fluxes (from which they originate) and that at- 
mospheric transport determines how much information on 
surface fluxes is represented within atmospheric CO2 con- 
centrations at any given location within the simulated atmo- 
sphere. To minimise the effects of short-term transport and to 
focus only on large seasonal variations, we conducted these 
analyses using monthly mean values. 

We investigated the representativeness of atmospheric 
CO2 concentrations simulated at TTA of surface CO2 flux. 
We compared the fraction of each ecosystem-specific net up- 
take CO2 tracer (e.g. for crop J, f crop ) simulated at TTA to the 
fraction of ecosystem-specific surface net CO2 uptake flux 
(e.g. for crop l /crop- where / indicates this is a flux). Each 
flux is the integral over the land surface included within the 
land surface mask detailed in Sect. 3.1. 
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Table 1. Tracer pools and definitions used by WRF-SPA. A non-interacting tracer does not have the potential to be exchanged with the land 
surface after its initial emission. An interacting tracer can be removed from the atmosphere as it represents a physical mass of CO 2 added to 
the atmosphere through respiration. 


Tracer 

Description 

Interacting tracer 

1 

Total CO 2 concentration; includes all sources and sinks of CO 2 for comparison to observations 

Yes 

2 

Forest net CO 2 uptake 

No 

3 

Anthropogenic emissions 

Yes 

4 

Forcings only, i.e. anthropogenic emissions, ocean sequestration, initial and lateral boundary conditions only 

No 

5 

Crop net CO 2 uptake 

No 

6 

Ocean sequestration 

No 

7 

Forest net CO 2 release 

Yes 

8 

Crop net CO 2 release 

Yes 

9 

Managed grassland net CO 2 release 

Yes 

10 

Other vegetation net CO 2 release 

Yes 

11 

Managed grassland net CO 2 uptake 

No 

12 

Other vegetation net CO 2 uptake 

No 


Y 4 Arop — 

Y 4 /crop 


. Arop 

4 Amp 4" 4 Aorest - ! - 4 Igrass - ! - 4 Ather 
4 /crop 

4 /cropd - 4 /forest - ! - 4 /grass - ! - 4 /other 


( 2 ) 

(3) 


Representativeness for any given ecosystem type is as- 
sumed to be when 


Y 4 Amp ^ Y 4 ./crop • (4) 

Moreover, this comparison provides an indication of 
whether the activity of a given ecosystem is over-represented 
(i.e. y 4 Arop > Y 4 /crop) or under-represented (i.e. y 4 
Arop < y 4 /crop) in simulated atmospheric CO 2 concentra- 
tions. Such information can help the interpretation of results 
from comparing ecosystem-specific CO 2 tracers and the sim- 
ulated atmospheric CO 2 concentrations. 

Investigation of seasonal variation is achieved through lin- 
ear regression analysis between the surface net CO 2 uptake 
flux and net uptake CO 2 tracer concentration for a given land 
cover type. As net uptake CO 2 tracers originate from the sim- 
ulated surface net CO 2 uptake flux, differences between sea- 
sonal variation of tracer concentrations and surface fluxes are 
due to how atmospheric transport relays variations in flux to 
TTA. To investigate interannual variation simulated to be de- 
tected at TTA, we assume that a change in surface net CO 2 
uptake flux for a given ecosystem should be reflected in the 
net uptake CO 2 tracers simulated at TTA. 


gradient from dominantly forested and peatland areas in the 
northwest to pasture in the central and southwest and arable 
cropland in the east. 

All meteorological data required, e.g. sea surface tem- 
perature (SST), soil initialisation, initial conditions and lat- 
eral boundary conditions, were taken from the Global Fore- 
casting System (GFS) reanalysis product (http://www.emc. 
ncep.noaa.gov/). GFS data are available at 1° x 1° longi- 
tude/latitude resolution with 6-hourly time steps (available 
from http://rda.ucar.edU/datasets/ds083.2/). The main fea- 
tures of the WRF model setup are presented in Table 2. 


3 Tall tower observations 

Observations of atmospheric CO 2 concentration are from 
TTA, a 222 m tower (observation height) near Dundee, Scot- 
land (56.56° N, 2.99° W). TTA is equipped for continu- 
ous measurement of atmospheric CO 2 concentrations, pro- 
ducing half-hourly observations which have been averaged 
to hourly timescales for comparison with WRF-SPA. TTA 
has been operational since the end of 2005 to the current 
date. Observations made at TTA have an accuracy limit of 
0. 1 ppm. TTA was part of the CHIOTTO network during 
the period of analysis reported here (EVK2-CT-2002-00163) 
and as such was fully integrated into the calibration and 
validation methodologies of that project. The data contin- 
ues to be quality controlled under the InGOS project (http: 
//www.ingos-infrastructure.eu/, accessed 9 December 2013, 
16:30 UTC). 


2.5 Model domain 


3.1 TTA footprint 


The WRF-SPA simulation is comprised of two grids in 
two-way nesting mode; the outer domain has a resolution 
of 18kmx 18 km and inner domain has a resolution of 
6 km x 6 km (Fig. 1). Scotland provides a highly complex 
topography and land use heterogeneity, with a longitudinal 


Currently there are no published assessment of TTA’s ob- 
servation footprint; however, the footprint of Mace Head, 
located on the west coast of Ireland, has been assessed in 
multiple studies (e.g. Henne et al., 2010; Rigby et al., 2011; 
Brunner et al., 2012). Mace Head is exposed to similar 
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Table 2. Parameter and model options used in WRF-SPA. 


Parameter 

WFA-SPA model options 

Basic equations 

Non-hydrostatic, compressible Advanced Research WRF (ARW) 

Radiative transfer scheme 

Rapid Radiative Transfer Model for GCMs (RRTMG) for both long- wave and short-wave 

Planetary boundary layer scheme 

Yonsei University 

Surface scheme 

Monin-Obukov 

Microphysics scheme 

WSM 3-class simple ice 

Cumulus parameterisation 

Grell 3-D ensemble scheme (coarse domain only) 

Nesting 

Two-way nesting 

Model time step 

Outer = 90 s, inner = 30 s 

Domain, resolution 

44 x47, 18 km 


48 x 54, 6 km 


35 vertical levels 

Domain centre 

56.63° N, 3.35° W 


meteorological conditions in northwest Europe, and there- 
fore we expect a similar footprint. Henne et al. (2010) calcu- 
lated a 12 h inversion, estimating the footprint of Mace Head 
to be the land surface within ~ 195 km of the tower. There- 
fore, if a similar footprint is assumed for TTA, ~ 98 % of the 
inner domain’s land surface is within the footprint of TTA. 
We use this estimate of footprint to mask the area of the land 
surface that is presented throughout this study. 

4 Results 

4.1 CO 2 time series at TTA 

We compared hourly observations of atmospheric CO 2 
concentrations from TTA (2006-2008) with both the 
WRF-SPA simulated total atmospheric CO 2 (R 2 — 0.67, 
rmse = 3.5ppm, bias = 0.58 ppm, linear regression) 
and “forcings only” CO 2 (R 2 = 0.71, rmse = 3.3ppm, 
bias = 0.82 ppm). These results suggest a slightly negative 
impact of including the modelled biosphere. However, the 
annual bias for total atmospheric CO 2 from WRF-SPA is 
lower than “forcings only” CO 2 , highlighting the impact 
of Scotland’s biosphere sink (Fig. 2a). Furthermore, the 
addition of biospheric fluxes captures diurnal variation in 
hourly observations, which is otherwise absent in “forcings 
only” CO 2 (Fig. 2ab). However, inclusion of biospheric 
fluxes results in an overestimation of night-time atmospheric 
CO 2 concentrations simulated at the tall tower (Fig. 2b). 
A comparison between atmospheric CO 2 concentrations 
observed at TTA and the CTE (2006-2007) and CT (2008) 
atmospheric inversion model (R 2 — 0.69, rmse = 3.5ppm, 
bias = 0.92 ppm) suggests comparable skill to WRF-SPA. 
However, CTE-CT does not simulate daily cycles in ob- 
served atmospheric CO 2 concentrations as well as WRF-SPA 
(Fig. 2b). As CTE-CT are available at 3-hourly intervals, the 
analysis was conducted against TTA observations averaged 
to a 3-hourly time step. 



01 June 2006 06 June 2006 11 June 2006 16 June 2006 21 June 2006 26 June 2006 


Fig. 2. Time series comparison between hourly observations of at- 
mospheric CO 2 concentrations made at TTA and WRF-SPA simu- 
lated total atmospheric CO 2 and “forcings only” CO 2 . Atmospheric 
CO 2 concentrations at 3-hourly time step from CTE-CT are also in- 
cluded in (b). (a) Shows that the simulated CO 2 time series (2006- 
2008) is mostly driven by forcings originating outside of the model 
domain, as indicated by “forcings only” CCF. (b) Shows an hourly 
(3 hourly for CTE-CT) time series for June 2006, highlighting that 
diurnal variation in simulated CO 2 is due to exchange with the 
biosphere within the simulated domain, as total atmospheric CO 2 
captures this variation. WRF-SPA modelled total atmospheric CO 2 
contains all model forcings and exchange with the simulated bio- 
sphere, while “forcings only” CO 2 does not include biospheric ex- 
change (i.e. total-biospheric fluxes comparison). 


The impact of the biosphere is more clearly seen at sea- 
sonal timescales using monthly means (Fig. 3). Total atmo- 
spheric CO 2 (R 2 — 0.96, rmse= 1.2 ppm, bias = 0.54 ppm), 
which includes biospheric exchange, shows improved sta- 
tistical agreement with observations compared to “forcings 
only” C0 2 (R 2 = 0.91, rmse=1.6ppm, bias = 0.71 ppm) 
and CTE-CT atmospheric CO 2 concentrations (R 2 — 0.94, 
rmse=1.5ppm, bias = 0.94 ppm). The monthly mean bias 
between total atmospheric CO 2 and observations is reduced 
for the majority of the comparison period relative to “forc- 
ings only” CO 2 and CTE-CT. The seasonal bias is reduced 
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Fig. 3. Time series of monthly mean residual (Model-Obs) between 
observed, CTE-CT and WRF-SPA simulated total atmospheric CO 2 
and “forcings only” concentrations. Highlights time periods during 
which the inclusion of the simulated biosphere results in a reduction 
in monthly mean bias. Error bars are ±1 standard error, accounting 
for temporal and spatial uncertainty only. 

in total atmospheric CO 2 by up to 2.8 ppm and 1.9 ppm be- 
tween March and June of each year compared to “forcings 
only” CO 2 and CTE-CT, respectively (Fig. 3). However, the 
modelled biosphere does not capture the observed seasonal 
minimum in atmospheric CO 2 concentrations, which occurs 
in July-August of each year (Figs. 2 and 3). During July- 
September total atmospheric CO 2 has a larger bias than both 
“forcings only” CO 2 and CTE-CT compared to observations. 
A larger positive bias in total atmospheric CO 2 than “forc- 
ings only” CO 2 indicates that modelled ecosystems within 
the footprint of the tall tower have become a net source of 
CO 2 at a time when they should remain a net sink (Fig. 3). 

4.2 Ecosystem contributions to atmospheric CO 2 
concentrations at TTA 

The dominant ecosystems simulated within the foot- 
print are forest, crop and managed grassland (Fig. 1). 
Over the validation period WRF-SPA simulated forest 
(— 2. 56 ± 0.05 tC ha -1 yr -1 , ± standard error, accounting 
for spatial and temporal uncertainty only) and man- 
aged grassland (— 0.48 ± 0.02 tC ha -1 yr -1 ) ecosystems to 
be mean annual sinks of carbon. Crop ecosystems 
(0.89 ± 0.01 tC ha -1 yr -1 ) were simulated to be a mean an- 
nual source of carbon. WRF-SPA estimates Scotland to 
be on average a carbon sink of — 0.99 ±0.04 tC ha -1 yr -1 . 
CTE-CT estimates Scotland to be a carbon source of 
+0.65 tCha -1 yr -1 . 

Net uptake and net release CO 2 tracers simulated at TTA 
suggest that cropland ecosystems have a distinct seasonal cy- 
cle from that of forests, managed grassland and “other” land 
cover types (Fig. 4). Managed grassland and “other” land 
covers are not included in the figure due to their contribu- 
tion to atmospheric CO 2 concentrations being small, never 
exceeding 0.7 ppm. Peak net uptake CO 2 tracer simulated 



Forest uptake 


01 Jan 2006 01 Jul 2006 01 Jan 2007 01 Jul 2007 01 Jan 2008 01 Jul 2008 

Date 

Fig. 4. Monthly mean mixing ratios for net uptake and net release 
CO 2 tracers for crop and forest ecosystems simulated at TTA. High- 
lights differences in detection of ecosystem processes at TTA, in 
particular the distinct seasonal cycle of cropland net uptake and 
net release CO 2 tracers compared to all other ecosystems. Man- 
aged grassland and “other” ecosystems are not included due to their 
small magnitude contributions, never exceeding 0.7 ppm. 

at TTA occurs 1 month earlier in crops than forest (except 
in 2006), while net release CO 2 tracer simulated at TTA 
for cropland shows a similar seasonality to all other ecosys- 
tems. Peak crop respiration coincides with crop harvest, a 
point in time when plant biomass has undergone senescence 
and has subsequently either been removed as part of har- 
vest processes or remains after harvest as residue added to 
the litter pool (mean simulated harvest day of year: 2006 
is 225 ± 17.9, 2007 is 229 ±21.4, 2008 is 229 ±21.5, stan- 
dard deviation accounting for spatial variability only). While 
crops contribute a similar amount of net release CO 2 tracer 
as forest during the growing season, during winter crop net 
release CO 2 tracer is less than half that of forest (Fig. 4). An- 
thropogenic CO 2 simulated at TTA is comparable in mag- 
nitude to that of CO 2 released by the biosphere. Also, an- 
thropogenic CO 2 does not display a strong seasonal trend 
(Fig. 4). 

Forest and crops dominate the net uptake CO 2 tracer sim- 
ulated at TTA (65-93 % of TTA tracer concentration and 72- 
91% of surface flux) (Fig. 5). On average, crop and man- 
aged grassland are over-represented in net uptake CO 2 trac- 
ers simulated at TTA by 3 % and 3.4%, respectively. Man- 
aged grassland represents on average just ~ 13 % of surface 
net CO 2 uptake flux compared to ~40 % for crops. Forest and 
“other” ecosystems are under-represented by 5 % and 1 .2 %, 
respectively. However, the over/under-representation varies 
at seasonal timescales, e.g. the largest under-representation 
of forests at TTA occurs between August 2006 and Jan- 
uary 2007 (21 %) while at other times atmospheric CO 2 sim- 
ulated at TTA is more representative. The bias towards crops 
is consistent with the spatial distribution of crops and for- 
est cover in relation to TTA’s location (Fig. 1). Crops is the 
dominant ecosystem, both in terms of net uptake CO 2 tracer 
and surface net CO 2 uptake flux during the growing season 
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Fig. 5. Comparison between monthly mean ecosystem-specific frac- 
tion of net uptake CO 2 tracers simulated at TTA and fraction of 
surface net CO 2 uptake flux. Were an ecosystems fraction of net 
uptake CO 2 tracer to be greater than the corresponding fraction of 
surface net CO 2 uptake flux, it would indicate that the ecosystem 
is over-represented in total atmospheric CO 2 concentrations. The 
reverse would indicate that the ecosystem was under-represented. 
Error bars are ± 1 standard error, accounting for temporal and spa- 
tial uncertainty only. 

(Fig. 5). After harvest (July of each year), forest becomes 
the dominant land cover for driving CO 2 exchange, as crop 
surface net CO 2 uptake flux declines due to senescence and 
removal of plants. 

4.3 Seasonal and interannual variation 

Net uptake CO 2 tracers simulated at TTA are able to explain 
the majority of seasonal variation in surface net CO 2 uptake 
for crops, forest, managed grassland and “other” land cov- 
ers (Table 3). The seasonal cycles in net uptake CO 2 trac- 
ers are more variable during the growing season (i.e. May- 
August), such that there is a mismatch between peak net up- 
take CO 2 tracers and surface net CO 2 uptake flux by ± one 
month (Fig. 6). Moreover, the amount of variation in surface 
net CO 2 flux explained by net uptake CO 2 tracers simulated 
at TTA varied between years (e.g. forest 2006 R 2 — 0.79 and 
2008 R 2 = 0.58). The rank order of net uptake CO 2 tracers 
simulated at TTA from each year does not correspond with 
the rank order of surface net CO 2 uptake flux for any ecosys- 
tem (Fig. 6). Interannual variation in mean annual surface 
net CO 2 uptake flux was ~ 9 % while interannual variation 
for mean annual net uptake CO 2 tracer at TTA was ~ 19 %. 

The annual prevailing wind direction over Scotland var- 
ied between years; in 2006 and 2008 the prevailing wind di- 
rection was broadly south/southwest, while in 2007 the pre- 
vailing wind direction was westerly. Moreover, the prevail- 
ing wind direction varied at seasonal timescales. During the 
peak growing season (May to August) there was consider- 
able variation (Fig. 7). In 2006 prevailing wind direction dur- 
ing the growing season varied between southerly and west- 
erly. In 2007 and 2008 there were periods of northerly and 


Table 3. Summary of R~ values from regression analysis of varia- 
tion in surface net CO 2 uptake flux explained by tall tower detected 
net CO 2 uptake tracers. A combined 2006 to 2008 period is pro- 
vided to give an indication of overall performance, while individual 
years allow for consideration of interannual variation in detection 
capability. 



Crop 

Forest 

Managed 

grassland 

“Other” 

2006-2008 

0.94 

0.72 

0.72 

0.77 

2006 

0.94 

0.76 

0.82 

0.85 

2007 

0.94 

0.74 

0.58 

0.71 

2008 

0.96 

0.58 

0.69 

0.70 


easterly winds, particularly during June; returning to more 
southwesterly directions by August in each year. This inter- 
annual and seasonal variation in wind direction would likely 
have impacted the detected footprint by TTA. 

5 Discussion 
5.1 CO 2 time series 

WRF-SPA demonstrated that it can recreate observations of 
atmospheric CO 2 concentrations (Fig. 2a and b). The dom- 
inant seasonal cycle reproduced by WRF-SPA is largely 
driven by forcings external to the modelled domain, i.e. the 
global signal from lateral boundary conditions (Fig. 2a), as 
indicated by the “forcings only” CO 2 tracer. Atmospheric 
CO 2 concentrations are underestimated during the winter; 
however, the bias is of a smaller magnitude in total atmo- 
spheric CO 2 than in “forcings only” CO 2 (Fig. 3). The un- 
derestimation during winter likely indicates that SPA un- 
derestimates net release of CO 2 flux (i.e. respiration) from 
the land surface. WRF-SPA has previously been validated 
against eddy covariance observations of NEE at forest, man- 
aged grassland and cropland sites, where forest and managed 
grassland NEE was overestimated during winter while crop- 
land was underestimated (Smallman et al., 2013). Given that 
net release CO 2 tracers simulated at TTA for crops is half the 
magnitude for forests, crops are a plausible candidate to ex- 
plain the wintertime underestimate in simulated atmospheric 
CO 2 concentrations (Fig. 4). Smallman et al. (2013) hypothe- 
sised that the underestimate in crop could be related to an un- 
derestimation of soil organic matter or the rate of soil organic 
matter turnover within the carbon model; however, there re- 
main several possibilities to be explored (e.g. ploughing). 

In contrast, CTE-CT continues to overestimate atmo- 
spheric CO 2 until late in the year (November/December), 
indicating that the inversion analysis continues to under- 
estimate Scotland’s carbon sink/overestimate carbon source 
during this period (Fig. 3). The WRF-SPA modelled bio- 
sphere generates diurnal cycles of realistic magnitude in the 
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Fig. 6. Seasonal and interannual comparison between monthly mean net uptake CO 2 tracer simulated at TTA for crop, forest, managed 
grassland and “other” (upper panels), and monthly sum surface net CO 2 uptake flux (lower panels). Note the different scales between 
ecosystem types. Error bars are ±1 standard error, accounting for temporal and spatial uncertainty only. 
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Fig. 7. Interannual comparison of growing season (May, June, July and August) prevailing wind direction at TTA. The wind rose shows the 
count of hourly wind directions simulated by WRF-SPA, where the direction indicated is the direction from which the wind is coming. 


modelled CO 2 time series (Fig. 2b) and reduces the seasonal 
bias seen in “forcing only” CO 2 (Fig. 3); however, nocturnal 
atmospheric CO 2 concentrations are overestimated. The noc- 
turnal overestimation of atmospheric CO 2 concentrations is 
likely due to an error in the YSU PBL scheme used in WRF- 
SPA. The error results in an overestimation of atmospheric 
eddy diffusivity under stable conditions, ultimately leading 
to a higher PBL (Hu et al., 2013). However, given that day- 
time CO 2 concentrations remain well simulated, it is unlikely 
that the nocturnal error persists into the well-mixed bound- 
ary layer due to rapid turnover of the atmosphere through 


nudging by LBCs. Moreover, WRF-SPA has been previously 
assessed against surface fluxes of heat, water and CO 2 , and 
daytime vertical profiles of atmospheric CO 2 concentrations 
where profile structure was well simulated, from which we 
can infer appropriate atmospheric transport (Smallman et al., 
2013). Furthermore, WRF-SPA’s performance is comparable 
to several studies that have compared observations of atmo- 
spheric CO 2 concentrations made at tall towers to high res- 
olution mesoscale model simulations (e.g. Ahmadov et al., 
2009; Tolk et al., 2009; Pillai et al., 2011). 
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Atmospheric CO 2 concentrations from CTE-CT were also 
compared to observation made at TTA, at both 3-hourly av- 
eraged (Fig. 2b) and monthly mean timescales (Fig. 3). The 
statistical comparison suggests little impact of high resolu- 
tion simulation using WRF-SPA, which contrasts with simi- 
lar comparisons between high and coarse horizontal resolu- 
tion models (e.g. Ahmadov et al., 2009). However, CTE-CT 
did not capture the observed diurnal cycle seen in TTA ob- 
servations as well as WRF-SPA (Fig. 2b). The mean bias be- 
tween observations of atmospheric CO 2 concentrations from 
TTA (which are not included in the CTE-CT atmospheric in- 
version) and CTE-CT is comparable to towers which were in- 
cluded in the inverse model (Peters et al., 2010). At seasonal 
timescales, CTE-CT tends to show a reduced bias compared 
with “forcings only” CO 2 ; however, the growing season bias 
remains larger than in total atmospheric CO 2 concentrations 
simulated by WRF-SPA (Fig. 3). 

The prevailing wind direction over the UK is southwest- 
erly, allowing the tower at Mace Head, Ireland, to provide an 
estimate of the background CO 2 concentration and to act as 
a boundary condition upwind of the air which passes over 
Scotland. Mace Head is used to provide a boundary con- 
dition in the CTE-CT atmospheric inversion (Peters et al., 
2010). Importantly, the bias between observations made at 
Mace Head and CTE-CT is small at +0.05 ppm (Peters et al., 
2010). Therefore, it can be inferred that the errors between 
modelled estimates of atmospheric CO 2 concentrations and 
observations is largely due to the simulation of surface ex- 
changes within the model domain presented here. 

5.2 Scotland’s carbon balance 

WRF-SPA’s estimate of carbon sink magnitude is ~5 fold 
greater than the official estimate of Scotland’s carbon sink by 
the UK National Atmospheric Emissions Inventory (NAEI), 
which estimates Scotland’s carbon balance, for CO 2 , to be 
-0.20 tC ha -1 yr -1 (Thomson et al., 2012). WRF-SPA does 
not account for a number of management impacts, such as 
biomass burning and land cover change. Excluding these 
fluxes, the NAEI estimate for Scotland’s carbon sink is 
— 0.44 tC ha -1 yr -1 . While WRF-SPA agrees with the NAEI 
that Scotland is a net sink of carbon, there appears to be a 
large discrepancy in the magnitude of the sink strength, the 
causes of which remain to be identified. However, as we cur- 
rently lack an error analysis of atmospheric CO 2 concentra- 
tions simulated by WRF-SPA, it remains unknown whether 
the discrepancy shown here is within errors. 

WRF-SPA simulated mean forest sequestration 
(— 2.56tCha _1 yr -1 ) is approximately double the estimates 
for UK wide (Cannell et al., 1999) and average European 
forest sequestration (Janssens et al., 2005; Luyssaert et al., 
2010). Scotland-specific estimates of forest sequestration are 
more similar to the simulations; Scotland-specific estimates 
range between ~ 1.8tCha _1 yr -1 (Thomson et al., 2012) 
and ~2.0tCha _1 yr -1 (Forestry Commission Scotland, 


2009). Forest activity is under-represented in atmospheric 
CO 2 concentrations simulated at TTA, as indicated by the 
lower fraction of net uptake CO 2 tracer simulated at TTA 
compared to the fraction of surface net CO 2 uptake flux 
originating from forest land cover. This under-representation 
of forest cover may explain why there is no apparent over- 
estimation of Scotland’s net carbon sink in the comparison 
between simulated CO 2 at TTA and observations. Grasslands 
were simulated to be a net carbon sink (—0.48 tC ha -1 yr -1 ) 
while croplands were simulated to be a net carbon source 
(0.89 tC ha -1 yr -1 ). Estimates of grassland carbon sink are 
more comparable with other estimates, which range between 
—0.69 tC hayr -1 (UK average, Janssens et al., 2005) and 
—0.15 tC hayr -1 (Scotland specific, Thomson et al., 2012). 
The WRF-SPA estimate of cropland source magnitude is 
also comparable with other UK wide (0.53 tC ha -1 yr -1 , 
Janssens et al., 2005) and Scotland-specific estimates 
(0.88 tC ha -1 yr -1 , Thomson et al., 2012). 

The simulated representation of arable cropland within 
WRF-SPA is likely to be responsible for the increase in the 
monthly mean bias between July-September in total atmo- 
spheric CO 2 (Fig. 3). Cropland net uptake CO 2 tracer simu- 
lated at TTA declines in magnitude concurrently with the in- 
crease in total atmospheric CO 2 bias in July (Figs. 3 and 4). 
In addition, the total atmospheric CO 2 bias exceeds the “forc- 
ings only” CO 2 bias in August as cropland respiration in- 
creases due to the input of litter from harvest. Above-ground 
carbon is removed as part of the harvest, leaving a fraction 
of above-ground carbon as surface residue and root carbon 
within the soil. Both the surface residue and root carbon are 
added to the litter carbon pool, which begin to decompose 
significantly increasing respiration from cropland. 

In Scotland, on average ~ 36 % of agricultural land is un- 
cultivated, including woodland patches, hedgerows and fal- 
low land (Scottish Government, 2012). However, WRF-SPA 
does not simulate uncultivated land associated with agricul- 
ture. These nonmodelled vegetative components are likely 
to be perennial systems, lack intensive management and as 
a result have a longer growing season. For example, forest 
and managed grassland ecosystems continue to have a sig- 
nificant level of surface CO 2 uptake flux for several months 
after cropland harvest (Fig. 6). Therefore, uncultivated sys- 
tems represent a significant contribution to the agricultural 
carbon balance at regional scales (Smith, 2004). Further de- 
velopment in the representation of agricultural land within 
LSMs is needed. For example, modelling at high spatial res- 
olutions may allow land cover maps to resolve some of this 
heterogeneity; alternatively, a tiling system could be used to 
represent this sub-grid heterogeneity. 

5.3 Representativeness and seasonal variation of TTA 
observations 

Cropland is most often the dominant ecosystem-specific 
net uptake CO 2 tracer simulated at TTA and the dominant 
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surface net CO 2 uptake flux (Fig. 6). Cropland is also frac- 
tionally over-represented at TTA compared to its surface net 
CO 2 uptake flux (Fig. 5). Over-representation of crops is ex- 
pected given the spatial distribution of the land cover types in 
relation to TTA (Fig. 1). Our results are consistent with other 
findings of both modelling and observational studies in this 
regard (Gerbig et al., 2009; Vermeulen et al., 2011; Lauvaux 
et al., 2012; Miles et al., 2012). Forests dominate the frac- 
tional activity after cropland harvest due to their continuing 
biological activity in the late summer (Figs. 5 and 6). 

Atmospheric CO 2 concentrations simulated at TTA con- 
tains significant seasonal information on ecosystems that are 
not adjacent to the tower (i.e. forest, managed grassland and 
“other”). Forest dominance of the fraction of net uptake CO 2 
tracer (Fig. 5) at TTA coincides with crop senescence and 
harvest (Fig. 4). The relatively small mismatches shown here 
seem likely to be explained by seasonal variation in tower 
footprint, as indicated by seasonal and interannual variation 
in prevailing wind direction (Fig. 7). Cropland is best rep- 
resented by net uptake CO 2 flux tracers (Table 3), which is 
expected given the local dominance already discussed. 

5.4 Interannual variation 

Interannual variation of the simulated seasonal cycles in sur- 
face net CO 2 uptake flux is poorly represented by net up- 
take CO 2 tracers simulated at TTA (Fig. 6). Interannual vari- 
ation of net uptake CO 2 tracers simulated at TTA is greater 
than interannual variation in modelled land surface net CO 2 
uptake flux. This suggests that interannual variation in at- 
mospheric transport due to year-to-year variation in weather, 
not variation in land surface net CO 2 uptake, is the dominant 
driver of interannual variation in tall tower observations for 
the years simulated here. This inference is supported by the 
interannual variation in wind direction, for example during 
the growing season (May, June, July and August) in 2008 
where there is a larger incidence of easterly winds altering 
the observation footprint of TTA (Fig. 7). This highlights the 
need for careful attention to atmospheric transport uncertain- 
ties and errors when carrying out atmospheric inversions. To 
detect a change in land surface activity, the magnitude of the 
change must be greater than the variation in detection due to 
transport. Alternatively, an extended network of tall towers is 
required to gain spatially explicit information on land surface 
exchange (Lauvaux et al., 2012). 

Net uptake CO 2 tracer concentrations for managed grass- 
land and “other” land covers simulated at TTA are less than 
0.7 ppm and 0.2 ppm, respectively. The accuracy limit for 
CO 2 detection of the equipment installed at TTA is 0. 1 ppm. 
This suggests that limited real world information is present 
in TTA observations for managed grassland and “other” land 
covers in the MODIS map. Therefore, it is also likely that 
TTA provides limited information in the real world for any 
ecosystem with limited activity or small spatial extent. In- 
effective detection of “other” vegetation is significant as 


“other” land cover types include Scotland’s uplands and peat- 
land areas. Upland and peatland areas are highly important 
given the significant amount of carbon stored as soil organic 
matter in these soils, estimated to contain >200tCha _1 
(Bradley et al., 2005). However, it should be noted that the 
MODIS land cover map used in WRF-SPA does contain 
errors; the upland and peatland cover for Scotland in the 
MODIS map is ~ 3 % while more detailed mapping efforts 
of Scotland have estimated uplands and peatlands to cover a 
larger area (~ 17 % Macaulay Land Use Research Institute, 
1993). 

5.5 Future work 

WRF-SPA estimates for ecosystem-specific mean annual se- 
questration are broadly reasonable; however, they should be 
considered with caution. This study does not estimate the 
SPA parameter uncertainties or uncertainties associated with 
atmospheric transport that may have a significant impact on 
the interpretation of sequestration estimates given here. As 
a result, estimates of ecosystem mean annual carbon seques- 
tration should be considered only as indicators of consistency 
with other estimates. Therefore, future work should involve 
an appropriate data-driven uncertainty analysis of SPA pa- 
rameters (i.e. data assimilation) and an attempt to assess at- 
mospheric transport uncertainties. 

WRF-SPA does not currently include a representation of 
forest management. Forest ecosystems are initialised with 
identical conditions that have been “spun up” into steady 
state. As a result, important differences in forest sequestra- 
tion due to age class distribution and lateral transport of 
carbon due to forest harvest are not included. In future, a 
more detailed representation of forest processes should be 
included. 

It remains to be investigated whether policy-relevant land 
cover management can be detected at TTA (e.g. afforesta- 
tion). WRF-SPA simulations presented here indicate that ob- 
servations made at TTA are unable to reliably detect interan- 
nual variation of ecosystems. However, tall towers are ex- 
pected to be used for monitoring the effects of land sur- 
face management aimed at mitigating climate change (ICOS, 
2012). Current Scottish government policy is to increase 
Scotland’s forest cover by 650 000 ha by 2050 (Forestry 
Commission Scotland, 2009). Through WRF-SPA the capa- 
bility of current observations to detect changes in Scotland’s 
regional carbon balance should be investigated. 

6 Conclusions 

Three specific questions were asked of WRF-SPA to inves- 
tigate atmospheric observations of CO 2 made mostly within 
the PBL from TTA, Scotland. 
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(i) Does WRF-SPA more accurately simulated observed 
atmospheric CO2 concentrations compared to a coarse reso- 
lution global atmospheric inversion model? WRF-SPA does 
more accurately simulate observed atmospheric CO2 con- 
centrations at TTA compared to CTE-CT. WRF-SPA bet- 
ter represents diurnal variation and a reduced bias between 
simulated atmospheric CO2 concentrations and observations, 
particularly during the growing season, (ii) Can ecosystem- 
specific CO2 tracers be used to inform on which ecosystem 
processes and land covers are responsible for observed vari- 
ations in atmospheric CO2 concentrations? Ecosystem spe- 
cific tracers have been successfully used to infer crops as 
being responsible for a increase in the bias between WRF- 
SPA simulated atmospheric CO2 concentrations at observa- 
tions post-harvest each year. Furthermore, we have hypoth- 
esised that the cause of the error is the lack of a representa- 
tion of uncultivated components of agricultural land not cur- 
rently parameterised for in WRF-SPA. (iii) Can observations 
made at TTA detect variation in ecosystem carbon uptake, 
for ecosystems within the footprint of TTA, at seasonal and 
interannual timescales? A majority of seasonal variation in 
surface net CO2 uptake flux is explained by net uptake CO2 
tracers for each ecosystem. However, the amount of varia- 
tion explained varied considerably between years. Moreover, 
interannual variation was not well captured, potentially due 
to seasonal and inter annual variation in the prevailing wind 
direction. However, for all other ecosystems interannual vari- 
ation in atmospheric transport due to year-to-year variation in 
weather had a greater impact on tall tower observations than 
interannual variation in surface uptake. 
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